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Abstract 

A lower bound is derived on the free energy (per unit volume) of a 
homogeneous Bose gas at density g and temperature T. In the dilute 
regime, i.e., when g <^ 1, where a denotes the scattering length of the 
pair-interaction potential, our bound differs to leading order from the 
expression for non-interacting particles by the term 47ra(2p^— [p— gc] + )- 
Here, gc{T) denotes the critical density for Bose-Einstein condensation 
(for the non- interacting gas), and [•]+ = max{ • ,0} denotes the posi- 
tive part. Our bound is uniform in the temperature up to temperatures 
of the order of the critical temperature, i.e., T ~ g^^^ or smaller. One 
of the key ingredients in the proof is the use of coherent states to 
extend the method introduced in [17] for estimating correlations to 
temperatures below the critical one. 



1 Introduction and Main Result 

The advance of experimental techniques for studying ultra-cold atomic gases 
has triggered numerous investigations on the properties of dilute quantum 
gases. From a mathematical point of view, several rigorous results have 
been obtained over the last few years. (See [8] for an overview.) The first of 
these, which has inspired much of the later work, was a study of the ground 
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state energy of a Bose gas with repulsive interaction at low density g. Per 
unit volume, it is given by 

eoig) = Airag^ + o{g^) for a^g < 1 (1.1) 

in three spatial dimensions. Here, a > denotes the scattering length of the 
interparticle interaction, and units are chosen such that h = 2m = 1, with 
m the mass of the particles. A lower bound on eo(^) of the correct form 
(1.1) was proved by Lieb and Yngvason in [11]. Much earlier, Dyson [2] had 
already proved an upper bound of the desired form, at least in the special 
case of hard-sphere particles. An extension of his calculation to arbitrary 
repulsive interaction potentials was given in [9]. 

The methods introduced in [11] have been extended to treat the case of 
fermions as well, for the study of both the ground state energy [7] and the 
free energy at positive temperature [16]. We are concerned here with the 
extension of (1.1) to positive temperature, at least as far as a lower bound is 
concerned. That is, our goal is to derive a lower bound on the free energy of a 
dilute Bose gas at density g and temperature T. Much of the complication in 
such an estimate is caused by the existence of a Bose-Einstein condensate for 
temperatures below some critical temperature. Although the existence of a 
condensate for interacting Bose gases has so far eluded a mathematical proof, 
its presence can easily be shown in the case of non-interacting particles. A 
short review of the Bose gas without interaction among the particles is given 
in Subsection 1.2 below. 

One of the main ingredients in our estimate is a method to quantify 
correlations present in the state of the interacting system. This method has 
been introduced in [17]; it does not immediately apply below the critical 
temperature for Bose-Einstcin condensation, however. We have been able 
to overcome this difficulty with the aid of coherent states. 

1.1 Definition of the Model 

We consider a system of bosons, confined to a three-dimensional flat 
torus of side lengths L, which we denote by A. The one-particle state space 
is thus iv^(A, dx), and the Hilbert space for the system is the symmetric N- 
fold tensor product TLn = ^sym(^^' d^ x), i.e., the space of square integrable 
functions of N variables that are invariant under exchange of any pair of 
variables. The Hamiltonian is given as 

N 

HN = J2-^i+ Hd{xi,Xj)). (1.2) 

1=1 l<i<j<N 
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Here, A denotes the Laplacian on A, and d{x,y) denotes the distance be- 
tween points X and y on the torus A. The particle interaction potential 

V : M-i- 1-^ M_(- U {00} is assumed to be a non-negative and measurable 
function. It is allowed to take the value +00 on a set of positive mea- 
sure, corresponding to hard sphere particles. In this case, the domain of 
the Hamiltonian has to be suitably restricted to functions that vanish on 
the set where the interaction potential is infinite. We assume that v has 
a finite range Rq, i.e., v(r) = for r > Rq. In particular, it has a finite 
scattering length, which we denote by a. We will recall the definition of a in 
Subsection 1.3 below. 

We note that in a concrete realization of A as the set [0, L]^ C M^, A is 
the Laplacian on [0, L]^ with periodic boundary conditions. Moreover, the 
distance d{x,y) is given as d{x, y) = min^-g^s \x — y — kL\. Note also that 
v{d{x, y)) = EfeGZS ^(1^ -y-kL\)\{L> 2Rq. 

The free energy (per unit volume) of the system at inverse temperature 
/3 = 1/T > and density ^ > is given by 

q) = ~ hm InTV^^ exp (-/3i?iv) , (1.3) 

where lim stands for the usual thermodynamic limit L ^ 00, N ^ 00 with 
g = iV/|A| fixed. Here, we denote the volume of A by |A| = L^. Existence of 
the thermodynamic limit in (1.3) can be shown by standard methods, see, 
e.g., [14, 15]. 

We arc interested in a bound on / in the case of a dilute gas, meaning 
that a^Q is small. The dimensionless parameter Pg^/^ is of order one (or 
larger), however. Note that sometimes in the literature the case of small g, 
but fixed a and /3, is understood with the term "dilute" . This corresponds 
to a high-temperature (classical) limit and is not what we want to study 
here. 

1.2 Ideal Bose Gas 

In the case of vanishing interaction potential {v = 0), the free energy can be 
evaluated explicitly. It is given as 

hW, a) = »up {„ + ^ * ,„ (1 - } . (1.4) 
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The supremum is uniquely attained at /x = /xo(/3, £>) = d/dgf{P, g) < 0. If ^ 
is bigger than the critical density 

the supremum is attained at = 0, whereas for g < gdP), it is attained at 
some Ho = fioiP, g) < 0. In particular, 

(2^ L e/^(p--^^°) - 1 = ^^^^^^ ■ ^^-^^ 

Note also that the scaling relation fo{P,g) = g^^^ foiPg'^^'^ A) holds for 
an ideal Bose gas. In particular, the dimensionless quantity /3^^/^ is the only 
relevant parameter. 

1.3 Scattering Length 

The scattering length of a potential v can be defined as follows (see Ap- 
pendix A in [12], or Appendix C in [8]): For R> Rq, 



Ana 
1-a/R 




dx {\Vct>{\x\)? + hv{\xmx\f) : 



(p:[0,R]<-^R+,(p{R) = lj. (1.7) 

For this definition to make sense, v need not necessarily be positive, one only 
has to assume that —A + (as an operator on L^(R^)) does not have any 
negative spectrum. We will restrict our attention to non-negative v, however. 
The infimum in (1.7) is attained uniquely. Moreover, the minimizer has a 
trivial dependence on R: for some function (py (independent of R) it can be 
written as (pi,{\x\)/(l)v{R)- Note that a is independent of R, and also that 
(/>„(|a;|) = 1 — a/\x\ for \x\ > Rq. 



1.4 Main Theorem 

Our main result is the following lower bound on the free energy, defined in 
(1.3). It gives a bound on the leading order correction, compared with a 
non-interacting gas, in the case of small a^g and fixed /3g^^^. 
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THEOREM 1 (Lower bound on free energy of dilute Bose gas). 

There is a function C : M+ ^ M+, uniformly bounded on bounded subsets of 
M+, and an a > such that 

fiP, q) > foiP, g) + 47ra (2^^ _ _ ^^(^)]2 ^ (i _ ^(i)) ^ (^ g) 

o(l)<C((/3^2/3)-l)(a^V3)a_ (19) 

Here, [•]+ = max{ • ,0} denotes the positive part. In the case of non- 
interacting particles, the expression [g — Qc{P)]+ is just the condensate den- 
sity. 

Remarks. 

1. Since C{t) is uniformly bounded for bounded t, our estimate is uniform 
in the parameter (Pg^^^)'^ as long as it stays bounded. I.e., our result 
is uniform in the temperature for temperatures not much greater than 
the critical temperature (for the non- interacting gas). In particular, 
we recover the result (1.1) in the zero temperature limit. The error 
term is worse, however; in [11], it was shown that the exponent a can 
be taken to be a = 3/17 at T = 0, whereas our proof shows that a can 
be chosen slightly larger than 0.00087 (independent of T). This value 
has no physical significance, however, it merely reflects the multitude 
of estimates needed to arrive at our result. 

2. The error term, o(l), in our lower bound depends on the interaction 
potential v, besides its scattering length a, only through its range 
Rq. This dependence could in principle be displayed explicitly. By 
cutting off the potential in a suitable way, one can then extend the 
result to infinite range potentials (with finite scattering length). See 
Appendix B in [9] for details. 

3. For Q < Qc{0) (i-e., above the critical temperature), the leading order 
correction term is given by Sirag^, compared with Airag^ at zero tem- 
perature. The additional factor 2 is an exchange effect; heuristically 
speaking, it is a result of the symmetrization of the wave functions. 
This symmetrization only plays a role if the particles are in different 
one-particle states, which they are essentially always above the critical 
temperature. Below the critical temperature, however, a macroscopic 
number of particles occupies the zero-momentum state; there is no 
exchange effect among these particles, which explains the subtraction 
of the square of the condensate density in (1.8). 
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4. We note that fo{(3, g) has a discontinuous third derivative with respect 
to Q at g = gdP) or, equivalently, a discontinuous third derivative with 
respect to T = at the critical temperature. Since the specific heat 
cv{P,g) can be expressed in terms of the free energy as cv{P,g) = 
—T(P/{dTy f{l3, g), it has a discontinuous derivative (with respect to 
r) at the critical temperature. The first order correction term in (1.8) 
has a discontinuous second derivative at this value. Considering only 
this term and neglecting higher order corrections, this would mean that 
the specific heat is actually discontinuous at the critical temperature. 

5. Our method applies also to particles with internal degrees of freedom, 
e.g., to particles with nonzero spin. For simplicity, we treat only the 
case of spinless particles here. 

6. Although we provide only a lower bound in this paper, one can expect 
that the second term in (1.8) gives the correct leading order correc- 
tion to the free energy (see, e.g., [4, Chapter 12.4]). To prove this, 
one has to derive an appropriate upper bound on /(/?, g), which has 
not yet been achieved, however. We note that a naive upper bound 
using first order perturbation theory yields (1.8) with Ana replaced by 
^ / dxv{\x\), which need not be finite, however (and is always strictly 
greater than Ana). 

The remainder of this paper is devoted to the proof of Theorem 1. The 
proof is quite lengthy and is split into several subsections. To guide the 
reader, we start every subsection with a short summary of what will be 
accomplished. 

Acknowledgments. It is a pleasure to thank Elliott Lieb and Jan Philip 
Solovej for many inspiring discussions. 

2 Proof of Theorem 1 

In the following, we find it convenient to think of A as the set [0, L]^ embed- 
ded in M^. We will also assume L to be large. In particular, L > 2Rq, but 
L will also be assumed to be large compared with several other parameters 
(with are independent of L) appearing below. This is justified since we are 
only interested in quantities in the thermodynamic limit L — oo. 

In many places in our proof, the Heaviside step function 9 will appear. 
We point out that we use the convention that 9 equals 1 at the origin, i.e., 
9{t) = for t < 0, and 9{t) = 1 for i > 0. 
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2.1 Reduction to Integrable Potentials 

Recall that we do not want to restrict our attention to interaction potentials 
that are integrable. For the Fock space treatment in the next subsection, 
it will be necessary that ?; has finite Fourier coefficients, however. As a 
first step, we will therefore replace the interaction potential f by a smaller 
potential v whose integral is bounded by some number Stti^. The scattering 
length of the new potential will be smaller than a, however. In the following 
lemma, we show that as long as ip is much greater than a, the change in the 
scattering length remains small. 

Lemma 1. Let v : R-|_ R-|_ U {00} have finite scattering length a. For 
any e > 0, there exists a v, with < v{r) < v{r) for all r, such that 
Jq°° dr r'^v{r) < 2ip, and such that the scattering length of v, denoted by a, 
satisfies 

a>a(l-^/a/^^{l-£) . (2.1.1) 

Proof. Without loss of generality, we may assume that ip > a. Let R = 
inf{s : dr r'^v{r) < 00}. We note that R is finite; in fact R < a. This 
follows from the fact that 2a > dr r'^v{r)\(/)y(r)\'^ , where (j)^ denotes the 
minimizer of (1.7) (for R = 00), as introduced in Subsection 1.3. Since it 
satisfies ^^(r) > 1 — a/r (see Appendix B in [12]), J°° drr^v{r) is finite for 
s > a. 

Assume first that drr'^v{r) > 2ip. The function s ^ drr'^v{r) is 
continuous for s > R. Wc can thus choose s > R such that drr'^v(r) = 
2(/9, and v{r) = v{r)6{r — s). 

To obtain an upper bound on a, we can use a trial function 0(r) = 
{(l}:^{r) — (j)^{s)s/r)9{r — s) in the variational principle (1.7). We note that (j) is 
a non-negative function, since <i)^{r) is monotone increasing in r [12]. By par- 
tial integration, using the variational equation — A^5j(|a;|) + |i;(|a;|)(/):5(|a;|) = 
0, we have 

47ra</ dx{\Vm)? + H\x\m\x\)\^) (2.1.2) 

drr^v{r)- (^</>g(s)- - 0^5(r)j . 

The last term is negative and can be dropped for an upper bound. To obtain 
an upper bound on s^(s), we note that (f>^{s) > 1 — a/s, and hence 
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For an upper bound on we use again the monotonicity of (f)v{r), which 

allows us to estimate 

drr'^v{r)(t)^{rf > (pvisfip. (2-1-4) 

This yields (f>^{s) < ^JaJIp. 

Altogether, we have thus shown that 

a < a + s(^{s) < a ( 1 + ) . (2.1.5) 

This proves (2.1.1) (with £ = 0) under the assumption that drr'^v{r) > 
2(p. 

Consider now the case when drr'^v{r) = 2ip — T for some T > 0. 
If R = 0, we can take v = v, and there is nothing to prove. Hence we 
can assume that R > 0. By definition, we have that drr'^v{r) = oo 

for any e > 0. Hence there exists a r (depending on T and e) such that 
Ir^i-c) '^^ min{t!(r), r} = T. We can then take 

v{r) for r > R 

v{r) = < min{u(r),r} for (1 - e)_R < r < (2.1.6) 
otherwise . 

Applying the same argument as in (2.1.2), with s = R, we have a < a + 
R(py{R). Now R < a, and (py{R{l — £)) < yfaff using the same argument as 
in (2.1.4), noting that v(r) = for r < R{\—e). Moreover, since | V</):^(|a;|)| < 
a/|xp, as shown in [9, Eq. (3.33)], |0^;(i?(l - e)) - <h{^)\ < £aR-^/{l - e), 
and thus 

a<a-^ + a.l^ . (2.1.7) 
1 — e y if 

This finishes the proof of the lemma. ■ 

As an example, consider the case of a pure hard sphere interaction, 
i.e., v{r) = oo for r < a, and v(r) = for r > a. In this case, we can 
choose v{r) = 6ipa^'^9(a — r). The scattering length of v is given by a = 
o(l — yja/ {'o^) tanh ■\jQip/a) in this case. Note that tanht < 1 for all t. In 
particular, a > a(l — ^/a/(6(p)). 

For a lower bound, we can simply replace v by v, i.e, we have > H^, 
with 

AT 

HN = J2-^i+ J2 ^d{xi,Xj)). (2.1.8) 

1=1 l<i<j<N 
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If we choose e < \fa[i^^ the error in the scattering length a is of the order 
y/a/ip. We will choose ip:^ a below. 

2.2 Fock Space 

In the following, it will be convenient to give up the restriction on the 
particle number and work in Fock space instead. This has the advantage 
that the problem of condensation can be dealt with with the aid of coherent 
states, which will be introduced in the next subsection. Our treatment differs 
slightly from the usual grand canonical ensemble since we do not simply 
introduce a chemical potential as Lagrange multiplier to control the number 
of particles, but we add a quadratic expression in N to the Hamiltonian. 
This gives a bitter control on the particle number. 

Let /Uq < be the chemical potential of the ideal Bose gas, which is the 
quantity that maximizes the expression in (1.4). Let = ^j^Hn be the 
bosonic Fock space over L^{A). Let aj, and Op denote the usual creation and 
annihilation operators of plane waves in A with wave functions L~^/^e~^P^. 
We define a Hamiltonian H on Fock space as 

M = T + V + K + /xoAr, (2.2.1) 

with 

T = Y.{p'-l^o)ala^ , V=^^^y(p)a^^a;_^a,a,, (2.2.2) 

P ' ' p,k,l 

and 

C 

K = 47ra— (N - . (2.2.3) 

Here and in the following, all sums are over p G ^Z^. The Fourier transform 
oiv is denoted by v, i.e., v{p) = Jj^dxv{\x\)e~'^^^ . It is uniformly bounded; 
in fact, \v{p)\ < v{Q) < Sinp, where p was introduced in the previous sub- 
section. The number operator J2p A>'^p denoted by N, whereas N is just 
a parameter. The parameter C is positive and will be chosen later on. 

The Hamiltonian H commutes with the number operator N, and can be 
thought of a direct sum of its restrictions to definite particle number. Note 
that the restriction to particle number N is just Hn, i-e, H = Hj^ on the 
sector of particle number N. This implies, in particular, that 

Tr^^ exp ( - (3Hn) < Tr^ exp (-/3M) . (2.2.4) 

We will proceed deriving an upper bound on the latter expression. 
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2.3 Coherent States 



To obtain an upper bound on the partition function TVjr exp(— /^H), we use 
the method of coherent states [10]. Effectively, this replaces the operators 
Op and Up by numbers. This can be viewed as a rigorous version of part of 
the Bogoliubov approximation, where one replaces the operators and Og 
by numbers. Such a replacement is particularly useful if the zero-mode is 
"macroscopically occupied", i.e, if aJaQ ~ |A|. We will use this method not 
only for p = 0, however, but for a whole range of momenta \p\ < Pc for some 
Pc ^ 0. Although not macroscopic, their occupation will be large enough to 
require this separate treatment. 

To be more precise, let us pick some Pc > and write !F = T< ® !F>, 
where J^^ and jr> denote the Fock spaces corresponding to the modes \p\ < 
Pc and \p\ > pc, respectively. Let M denote the number of p G with 
\p\ < Pc- As shown in [10], the Berezin-Lieb inequality [1, 5] implies that 

TV^ exp(-/3M) < / d^z Tr^> exp ( - pm^{z)) . (2.3.1) 

Here, z denotes the vector (zi, . . . , zm) G C^, d^z = YifLi and dz = 
7r~^dxdy with x = ^{z), y = Moreover, EI'^(z) is the upper symbol 

of the operator H. It is an operator on jr>, parametrized by z, and can be 
written in the following way. Let \z) G J^<: denote the coherent state 

\z) = exp (E|p|<p,%4 - 4^p) |0) ^ ' (2-3-2) 

with |0) the vacuum in the Fock space J^<. Then the lower symbol of H is 
given by M.s{z) = {z\M.\z). Since ap\z) = Zp\z), the lower symbol is obtained 
from the expression (2.2.1) by simply replacing all the ap by Zp and the aj, 
by z* for all \p\ < Pc- The upper symbol can be obtained from the lower 
symbol by replacing jzpj^ by \zp\'^ — 1, for instance, and similarly with other 
polynomials in Zp-, see, e.g., [10] for details. We can then write IP{z) in the 

following way. Denoting by Ns(-?) = \z\'^ + Yl\p\>pc ^p^p lower symbol of 
the number operator, we have 

m^{z) = m^iz) - AM(z) , (2.3.3) 
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with 

^ 7(0) {2MNs{z) - M^) 



|p|<Pc ' ' 



+ 2 Yl v{l-k)ala,+ ^ vil - k) {2\zkf - l) 

\l\<PcAk\>Pc \l\<pc,\k\<pc 

+ [2|^2 ^ M(2Ns(z) - 2Ar - M)] . (2.3.4) 

Here, we have used that v{p) = v{—p). 

Since ■u is a non-negative function, \v{p)\ < v{0) < Sir^p for aU p. Hence 
we obtain the bound 

Am{z) < M{pI - ^o) + ^^N,(z) + ^ [|.f + M(N,(z) - AT)] . 

(2.3.5) 

(Here, we have used again the positivity of v.) Denoting by ]Ks(z) = (z|]fC|z) 
the lower symbol of IC (and, similarly, for T and V below), we have 

47raC ,^ ,^^9 . |-,|2\ ^ 47raC 



^.{^ = ((Ns(^) - NY + \z\') > ^ mz) - NY . (2.3.6) 

Wc can use part of IKs(i') to estimate —AH from below independently of z. 
More precisely, we have 

^Ks(z) - AM{z) 

> _M(p? - ,o) - ^(2.M + HO - 32™C<^ (l + ' 

= -Z(^) . (2.3.7) 

Note that M ~ Pcl-^l thermodynamic limit. We will choose the 

parameters pc, (f and C such that Z^^^ <C |A|a^^ for small ag^^^. 
With the definition 

Fg{p) = ln'IV^> exp ( - /3(T,(z) + W,{z) + ^K,{z))) , (2.3.8) 
(2.3.1) and the estimates above imply that 
-^InTV^ exp(-/3H) > i^oN-^ln I d^z exp (-/3F^(/3)) -Z^^^ . (2.3.9) 
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Hence it remains to derive a lower bound on -FV(/?) . 

Let r^- denote the Gibbs state of Ts(z)+Vs(^ + ^lfCs(£) on for inverse 
temperature /3. Let Hq = |0)(0| denote the vacuum state in Denoting 
by the state = U{z)IloU{z)^ ® on the full Fock space J^, we can 
write 

=Tr^[(T + V+^K)T^ - ^S(T^) . (2.3.10) 
Here, ^(F) = — Trj^FlnF denotes the von-Neumann entropy. 

2.4 Relative Entropy and A Priori Bounds 

In the following, we want to derive a lower bound on Fg{(]). Although we 
do not have an upper bound available, we can assume an appropriate upper 
bound without loss of generality; if the assumption is not satisfied, there 
is nothing to prove (as far as a lower bound in concerned). This upper 
bound can be formulated as a bound on the relative entropy between the 
state = U{z)IloU{z)^ (E> F^ defined above and a simple reference state 
(describing non-interacting particles). Together with a bound on the total 
number of particles, this estimate on the relative entropy contains all the 
information we need in order to prove the desired properties of the state 
that will allow us to derive a lower bound on (2.3.10). 

Wc note that, for any state F of the form F = U{z)IloU{z)^ (8) F^ for 
some state F^ on .F> , 

IV^[TF] - ^S(T) = TV^, [T,(z)F>] - ^S{T>) 

> ln'IV^> exp ( - (3T,{z)) . (2.4.1) 

In fact, the difference between the right and left sides of (2.4.1) is given by 
P~^S{r,UQ), where S denotes the relative entropy. For two general states 
F and F' on Fock space, it is given by 

S'(F,F') = TrjrF(lnF-lnF') . (2.4.2) 

Note that the relative entropy is a non-negative quantity. The state JIq is 
given by ilf = U{z)UoU{z)^ F°, where F° is the Gibbs state of Ts{z) on 
T> (which is independent of z). 
For F = T^, we have 

5(T^ ) = Tr^> F,-( In F,- - In F°) = 5(F^, F°) . (2.4.3) 
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Prom these considerations, together with the positivity of V, we conclude 
that (2.3.10) is bounded from below by 

Fm > -;^lnTV^^ exp ( - ^T,{z)) + iTr^[KT^^ + ^S{Tg,T^) . (2.4.4) 

Hence we can distinguish the following two cases: 

A) The following lower bound on Fg{j3) holds: 

Fm > ~ lnTV^> exp ( - (3T,{z)) + 87r|A|S^2 _ (2.4.5) 

B) Inequality (2.4.5) is false, in which case 

S{Ts, r°) < 87r|A|S/?^2 (2.4.6) 

and 

TV^[KT^ < 167r|A|S^2_ (2.4.7) 

From now on, will consider case B, i.e., we will assume (2.4.6) and (2.4.7) 
to hold. The lower bound we will derive on Fg{P) below will actually be 
worse then the bound (2.4.5) above; i.e., the bound in case B holds in any 
case, irrespective of whether the assumptions (2.4.6) and (2.4.7) actually 
hold. 

Although the relative entropy does not define a metric, it measures the 
difference between two states in a certain sense. In particular, it dominates 
the trace norm [13, Thm. 1.15]: 

5(r,r') > 5l|r-r'||i. (2.4.8) 

This inequality is a special case of the fact that the relative entropy de- 
creases under completely positive trace-preserving (CPT) maps. In fact, 
inequality (2.4.8) can be obtained using monotonicity under the CPT map 
r I— > Tr_7r[Pr] © TrjF[(l — P)T], where P is the projection onto the subspace 
where P - P' > 0. 

Although we have the upper bound (2.4.6) on the relative entropy, in- 
equality (2.4.8) is of no use for us since the relative entropy is of the order 
of the volume of the system, while the right side of (2.4.8) never exceeds 2. 
To make use of (2.4.8), we must not look at the state on the full Pock space 
(over the whole volume) but rather on its restriction to a small subvolume. 
We do this in Subsection 2.8 below. Again, the monotonicity of the relative 
entropy will be used in an essential way. 
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We note that (2.4.7) implies the following simple upper bound on \z\^. 
Prom (2.2.3) and (2.4.7), 

|il2_iv<Tr^[(N-iV)T^^ < (TV^[(N-iV)2T^"]y^^ < -A|A|^, (2.4.9) 

and hence 

We will choose C 2> 1 below. 
2.5 Replacing Vacuum 

In the following, we want to derive a lower bound on the expectation value 
of the interaction energy V in the state T^, i.e., on 

TV^ [VT^^ = Tr^> [V,(f)r,-] . (2.5.1) 

For reasons that will be explained later (sec Subsection 2.13), we find it 
necessary to replace the vacuum ITo on in the definition of the state = 
U{z)IloU{z)^ ® by a more general quasi-free state. In this subsection, 
we show that such a replacement can be accomplished without significant 
errors. 

Let n denote a (particle- number conserving) quasi free state on J^^. It is 
completely determined by its one-particle density matrix, which we choose 
to be given as 

\p\<Pc 

Here, \p) G L'^{^) denotes a plane wave of momentum p. We denote the trace 
of TT by P = Yl\p\<pc ^p- '^TT denote the state = U {z)IlU{z)^ (g) Tg on 
T. We want to derive an upper bound on the difference 

TV^[V(T^ - T^^] = TV^[V(i7(£) (H - Hq) U{z)^ ® Tg)] . (2.5.3) 
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A simple calculation yields 

+ 21X1 ^ - I) [TTkTTl + 2\Zk\'^Tri] 

' ' \k\<Pc,\l\<Pc 

+ m Yl ~ '^''■^> 

' ' |fc|<Pc,|/|>Pc 

< ^ (^p2 ^ 2P Tr^ [NT^^ ) . (2.5.4) 

Here we have used again that \v{k)\ < v{0) < Sirip. It follows easily from 
(2.4.7) (compare with (2.4.9)) that Tr:r[NT^ < N{l + 2/VC). Hence we 
obtain from (2.5.4) that 



with 



TV^[VT^^ > Tr^[VT^] - Z^^) , (2.5.5) 

|A| |A| V VCJ 

RecaU that C > 1 and <y? > a. Hence Z^^) < \A\aQ^ as long as (pP < aN. 
Note that the effect of the replacement of by T| on the kinetic energy 



IS 



TV^[TT-^ = TV^[TT^] - ^ (p2 _ f,^)^^ . (2.5.7) 

\p\<Pc 

We have thus obtained the lower bound 

F,{(3) > TV^[(T + V)Tf] - ^S{T^) 

|p|<Pc 



2.6 Dyson Lemma 

Since the interaction potentials in V are very short range and strong (com- 
pared with the average kinetic energy per particle), we cannot directly obtain 
information on the expectation value of V in the state T^. In fact, we can- 
not even expect that it yields the desired correction to the free energy, since 
part of the interaction energy leading to the second term in (1.8) is actually 
kinetic energy! Hence we will first derive a lower bound on V in terms of 
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"softer" and longer ranged potentials, with the aid of part of the kinetic 
energy. More precisely, wc will use only the high momentum part of the 
kinetic energy for this task, since this is the relevant part contributing to 
the interaction energy. The appropriate lemma to achieve this was derived 
in [7] ; part of the idea for such an estimate is already contained in the paper 
by Dyson [2]. For this reason, we refer to this estimate as "Dyson Lemma". 

Our goal is to derive an appropriate lower bound on the Hamiltonian 
T + V. Let X : IK be a radial function, < x{p) < 1, and let 

^(^) = ^E(1-X(p))e"'^"- (2-6.1) 

We assume that x{p) — 1 as |p| — oo sufficiently fast such that h G L^{A.) D 
L°°{A). For some L/2> R> Ro, let 

fnix) = sup \h{x -y)- h{x)\ , (2.6.2) 
\y\<R 



and 



wr{x) = ^fR{x) [ dyfRiy) . (2.6.3) 
Ja 



Note that wr is a periodic function on M^, with period L. 

Let Ur : M+ i— R_,_ be a non-negative function that is supported in the 

interval [Rq,R], and satisfies dtt^UR{t) < 1. The following is a simple 
extension of Lemma 4 (and Corollary 1) in [7]. The proof follows closely the 
one in [7]. For completeness, we present it in the appendix. 

Lemma 2. Let yi, . . . ,y„ denote n points in A and, for x G A, let yt^j<i{x) 
denote the nearest neighbor of x among the points yj, i.e., the y^ minimizing 
d{x,yj) among all yj. We then have, for any e > 0, 



-Vx(p)^V + 5^i;(d(x,yi)) > {l-e)aUR{d{x,yt^f^{x)))-^-WR{x-yi). 
1=1 1=1 

(2.6.4) 

Here, the operator —Vxip)"^^ stands for Y^pp'^x(p)^\p) {p\- 

We note that yf^j<i{x) is well defined except on a set of measure zero. 
Compared with Lemma 4 in [7], the main differences are the boundary 

conditions used, and the fact that wc do not demand a minimal distance 
between the points y^. In [7], it was assumed that d{yi,yj) > 2R for i ^ j, 
in which case UR{d{x,yj^]^{x))) = J2i UR{d{x,yj)). Note also that only the 
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inequality / dtt^U{t) < 1 is needed for the estimate, not equality, as stated 
in [7]. 

We will use Lemma 2 for a lower bound to the operator T + V on .7^. 
Note that the restriction of this operator to the sector of n particles is just 
Hn, defined in (2.1.8). We write 



(2.6.5) 



and apply the estimate (2.6.4) to each term in square brackets, for fixed j 

and fixed positions of the Xj, i 7^ j. We want to keep a part of the kinetic 
energy for later use, however. To this end, we pick some < x < 1, and 
write 

p2 = _ (1 _ x)x(p)') + (1 - H)p'x{vf ■ (2.6.6) 

We split the kinetic energy in the Hamiltonian (2.6.5) accordingly, and apply 
(2.6.4) to the last part. Using also the positivity of the v, we thus obtain, 
for any subset Jj C {1, . . . — 1, j + 1, . . . , n}, 

>-V,(l-(l-x)x(p,)2)V,- ^ (2.6.7) 

e 



+ (1 - £)(1 - K)aUR{d{xj,x^^{xj))) - - X] WR{xj - Xi) . 



Here wc denoted by x^-j^{xj) the nearest neighbor of Xj among the points 

X J 1/ ^ J j • 

Our choice of Jj will depend on the positions Xj, i 7^ j. We want to 
choose it in such a way that d{xi,Xk) > R/5 if I G Jj and k G Jj. Moreover, 
we want the set to be maximal, in the sense that if Z Jj, then there exists 
a /c € Jj such that d{xi,Xk) < R/5. These properties of Jj will be used in 
an essential way in Subsections 2.9 and 2.10 below. 

There is no unique choice of Jj satisfying these criteria. One way to con- 
struct it is the following. We first pick all i corresponding to those Xj whose 
distance to the nearest neighbor (among all the other x^, k 7^ is greater 
or equal to R/5. Secondly, going through the list {xi, . . . , Xj_i, Xj+i, . . . , x„} 
one by one, we add i to the list if d{xi,Xj) > R/5 for all j already in the 
list. This last procedure depends on the ordering of the Xj, and hence the 
resulting Jj will depend on this ordering. The resulting interaction poten- 
tial in (2.6.7) will thus be not symmetric in the particle coordinates. This is 
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of no importance, however, since we will take the expectation value of the 
resulting operator only in symmetric (bosonic) states anyway. 

The set Jj is chosen in order to satisfy the following requirements. On the 
one hand, we want the particles keep a certain minimal distance, -R/5; this 
is necessary in order to control the error terms coming from the potentials 
wr. We do not have sufficient control on the two-particle density to control 
these terms if all the particle configurations were taken into account. On the 
other hand, we want the balls of radius R centered at the particle coordinates 
to be able to overlap sufficiently much, such that the desired lower bound 
can be obtained. We note that we want to derive a lower bound which 
is independent of z] for certain values of z, however, the system may be 
far from being homogeneous and particles may cluster in a relatively small 
volume. We want to ensure that there is still sufficient interaction among 
them. 

2.7 Filling the Holes 

One defect of Lemma 2 above is that the resulting interaction potential Ur 
is supported outside a ball of radius Rq, which is the range of v. For our 
estimates in Subsection 2.9, it will be convenient to have a specific Ur which 
is, in particular, positive definite and hence should not have a "hole" at the 
origin. We will show in this subsection that one can easily add the missing 
part to Ur, at the expense of only a small amount of kinetic energy. 

We start with the description of our choice of Ur. Let j : R+ — M+ 
denote the "hat function" 

m = —[ dye{l-\y\)9{l-\y-et\) (2.7.1) 

for some unit vector e G M^. Note that j is supported in the interval [0, 1], 
and Jq dtt^j(t) = 1. An explicit computation yields 

j{t) = 12{t + 2)[l-t]l. (2.7.2) 

Our desired interaction potential will be UR{t) = R~^j{t/R). We will thus 
choose UR{t) = UR{t)e{t - Ro) in (2.6.4). 

In the following, it will be convenient to work with Ur instead of Ur. I.e., 
we would like the add the missing part Ur{ ■ )9{Ro— ■ ) to the interaction. In 
order to achieve this, we use the following lemma. It is an easy consequence 
of the definition of the scattering length, given in (1.7). 
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Lemma 3. Let yi,---,yn denote n points in k, with d{yi,yj) > R/5 for 
i^j. LetO<X< 7r/2, and Rq < R/10. Then 



" 1=1 



> - 



^^i3(^^-i)|:wio-.(x.„.)). 



Proof. It suffices to prove that 

i2 



J|a:|<R/10 -fl-O -'kl<-Ro 

for any function (j) G i/^. In fact, it is enough to prove (2.7.4) for radial 
functions. Note that the scattering length of the potential 2\^Rq^9{Ro — • ) 
is given by i?o(l ~ A~^tanA). Hence, for any Rq < s < R/IO, 

[ |V^(|x|)|^ - ^ / \H\x\)f > -4^Ro - l) \4>{s)f . 

J\x\<s -n-o J\x\<Ro \ ^ / 

(2.7.5) 

Eq. (2.7.4) follows by multiplying this inequality by and integrating s 
between Rq and R/10. ■ 

Let A = 7r/4 for concreteness. Recall that d{xi,Xk) > R/5 for i,k £ Jj. 
Since Unit) < j{Q)/R^ = 24/R^, and UR{t) > j{l/10)/R^ for t < R/10, 
this lemma implies, in particular, that 



{Un - UR)idixj,xi^^{xj))) < -^^-^Aj (2.7.6) 
+ (W " (i^/4°_^3,-(l/\o) ^^(^(^.--^NNfe))) • 



Define 



a' = a(l - e)(l -x)(l- (4 - tt) ——^ ^ , , ^ (2.7.7) 

^ '\ (7r/4)3^ ' {R/10f-Rlj{l/10)) ^ ' 

and 

24a (4Ro)^ 
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In the following, we will choose x S> aE^/R^ and hence, in particular, 
x' > 0. Combining the estimates (2.6.7) and (2.7.6) and applying them in 
each sector of particle number n, we obtain the inequality 



T + V > T'^ + ' 



(2.7.9) 



where 



T'^ = ^£(p)atap , e{p) = x'p^ + {l-x)p\l-x{p?)-fio, (2.7.10) 



and W is, in each sector with particle number ra, given by the (symmetriza- 
tion of the) multiplication operator 



n I- ~ 

^ a'UR{d{xj,Xj^^{xj))) ^ WR{xj - Xi) 



(2.7.11) 



Note again that the set Jj depends on all the particle coordinates Xi, i ^ j. 
We now describe our choice of the kinetic energy cutoff x- Let : i— 

M+ be a smooth radial function with i^(j)) = for \p\ < 1, iy{p) = 1 for 
\p\ > 2, and < vip) < 1 in-between. For some s > R we choose 



X{p) = v{sp) 



(2.7.12) 



We will choose < 1/s below. This imphes, in particular, that e(p), defined 
in (2.7.10) above, is equal to (1 — x+x')p^ — /xq for \p\ < pc- Hence (compare 
with (2.5.7)) 

IV^[rn] = Tr^[rT^^ + ^ ((1 - x + x')p' - /xo)% • (2.7.13) 

\p\<Pc 

Using the fact that 

TV^[T^T-^ - ^S{T^) > ~ InTV^^ exp ( - /3T,^(/)) , (2.7.14) 
we conclude from (2.5.8), (2.7.9) and (2.7.13) that 



Fm > --^lnIV^>exp(-/3T^(z)) + TV^[WTf] + 1tV^[KT^"] 



-(x-x') P%-Z^^^- 

\P\<PC 



(2.7.15) 
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Note that the first term on the right side of (2.7.15) can be computed ex- 
phcitly. It is given by 

--^ InTV^^ exp ( - /3T^(z)) = ((1 - x + V)p' - ^io) \zp\^ 

\P\<PC 

+ ^ 5; In (1 - exp ( - Pe{p)) . (2.7.16) 

|p|>Pc 

In the following, we will derive a lower bound on the expectation value of 
W in the state T£. 



2.8 Localization of Relative Entropy 

Our next task is to give a lower bound on Trjr [WT^] . For that matter, we 
will show that wc can replace the unknown state in = U {z)IlU{z)^^r^ 
by the quasi-free state T^, which is the Gibbs state for the kinetic energy 
Ts{z). The error in doing so will be controlled by the upper bound on the 
relative entropy, Eq. (2.4.6). In order to do this, we have to obtain a "local" 
version of this bound. 

Consider the quasi- free state = U(^r°. Its one particle density matrix 
is given by 

p 

with £{p) = (3{p'^ — fio) for \p\ > pc, and £{p) = ln(l -|- I/tt^) for \p\ < pc- 
Let 77 : H- > M be a function with the following properties: 

• 77(0) = 1, and r]{x) = for |x| > 1 

• V{p) = J dx r]{x)e~''P'' > for ah p G M^. 

An appropriate rj can, for instance, we obtained by convolving a smooth 
function supported in a ball of radius ^ with itself. Given such a function r), 
we define rib{x) = r]{x/h) for some h < L/2. Moreover, with a slight abuse 
of notation, we define a one-particle density matrix Uf, on H by the kernel 

i^bix, y) = to-nix, y)r]b{d{x, y)) . (2.8.2) 

Note that this defines a positive operator, with plane waves as cigcnstatos. 
Note also that \LOii{x,y)\ < \iOTr{x,y)\ since \r]i,\ < 1. We denote by Of, 
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the corresponding (particle number conserving) quasi-free state on T ^ and 

= U (z^Q.}^ [z)'^ . Let also denote = ijJ\X^^x) = uJt^(x,x) the one- 
particle density of fit (which is independent of x). Abusing the notation 
even more, we shall sometimes write u}b(x, y) = ijJh(x — y) if no confusion can 
arise. 

For r < L/2, let Xr,^{ ■ ) = 0{r—d{ ■ , ^)) denote the characteristic function 
of all ball of radius r centered at ^ G A. The function Xr,g defines a projection 
on the one-particle space TL = L?{h.), and hence the Fock space over Ji 
can be thought of as a tensor product of a Fock space over Xr,{H and a Fock 
space over the complement. States on T can thus be restricted to the Fock 
space over Xr,(,'H, simply be taking the partial trace over the other factor. 
We denote such a restriction of a state F by F^^ ^ . 

For d{^, C) > 2r, Xr,^ + XrX defines a projection on H. Note that since 
Ub{x,y) vanishes if d{x,y) > b, we have that 

^b,Xr,i+XrX = ^b,Xr,( ^ ^b,XrX (2.8.3) 

if d{^, C) > 2r -I- b. This follows simply from the fact that the one parti- 
cle density matrix of ^^6,^^^+^^^ is given by (xr,e + Xr,c)^b{Xr,^ + Xr,c) = 
Xr,$,'^bXr,( ~^ Xr.(^bXr.c,- The Same factorization property (2.8.3) is obviously 
true with $7^ replaced by = U{z)Qi,U{z)^ since the unitary U{z) has 
the same product structure. As in [17, Sect. 5.1], we have the following 
superadditivity property of the relative entropy. 

Lemma 4. Let Xi, < i < k, denote k mutually orthogonal projections on 
TC . Let Q be a state on T which factorizes under restrictions as . = 
^i^lxi- Then, for any state T, 

s{T,n)>J2s{rx„nx,). (2.8.4) 

i 

We note that the lemma applies, in particular, to a (particle number 
conserving) quasi-free state O whose one-particle density matrix u satisfies 
XiOjXj = for ? 7^ J. We emphasize that the factorization property of is 
crucial; in general, the relative entropy need not be superadditive. This is 
the reason for introducing the cutoff b. 

Proof. Let X denote the projection X = Yli-^i- The relative entropy de- 
creases under restrictions [6, 13], i.e., 

-s(F, n) > s{rx. nx) = s{^x,0^nx,) 

= ^ siFx^nx,) + Yl s{rx,) - sivx) . (2.8.5) 

i i 
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The last two terms together are positive because of subadditivity of the 
von-Neumann entropy. ■ 



We take the Xi to be the multiplication operators by characteristic func- 
tions of balls of radius r, separated a distance 2b. By averaging over the 
position of the balls, Lemma 4 implies that, for any b >2r such that L/{2b) 
is a positive integer, and for any state F, 

s{tM) > j^l^d^s{r^^^^,nl^j . (2.8.6) 

We apply this inequality to the state T = Tf = U{z)UU{zy (g) Tg. 

We remark that that restriction of L/(2b) being an integer will be of no 
concern to us, since we are interested in the thermodynamic limit L — > oo, 
with b independent of L. 

We can now apply inequality (2.4.8) to the right side of (2.8.6). Using 
the Schwarz inequality for the integration over ^, we thus obtain 

IITJ^^^ - nl^JU < 4 (b'lm^iM))'^' (2.8.7) 

for any r < 6/2. Note that ^(T^, $7^) = S(T.,^, fl^) since the relative entropy 
is invariant under unitary transformations. Were it not for the cutoff b, we 
could use (2.4.6) to bound the right side of (2.8.7). We will estimate the 
effect of the cutoff in Subsection 2.13. 

2.9 Interaction Energy, Part 1 

The next step is to derive a lower bound on Ttjf [WT^] . The main input 
will be the bound (2.8.7) derived in the previous subsection. We split the 
estimate into three parts. First, we give a lower bound on the expectation 
value of the terms containing Ur in (2.7.11). In the next subsection, we 
bound the remaining energy containing the terms wr. Finally, we combine 
the two estimates in Subsection 2.11. One of the difficulties in our estimates 
results form the fact that z is rather arbitrary, and hence the system can be 
far from being homogeneous. 
From (2.7.11), we can write 

W = Wi - Wa , (2.9.1) 

where 

oo n 

Wi = j;a't7K(d(x,-,4^(x,))) (2.9.2) 

n=0 j=l 



23 



and 

oo n ~ 

= © E E - • (2.9.3) 

n=0 j=l ieJj 

We start by giving a lower bound on the expectation value of Wi in the state 
T^. Recall that is defined after Eq. (2.5.2) as = U{z)UU{z)^ (g) Tg. 
According to the decomposition (2.7.1), we can write Ur as 

URidix, y)) = ^ [ e{R/2 - di^, x))e{R/2 - d{^, y)) . (2.9.4) 

This gives rise to a corresponding decomposition of Wi, which we write as 

144a' 



ttR^ 



[ d^wiO. (2.9.5) 
J A 



For r > 0, let n,.^^ denote the operator that counts the number of particles 
inside a ball of radius r centered at ^ G A. It is the second quantization of 
the multiplication operator Xr,^{ ■ ) = G{r — d{^, ■ )) on -C'^(A). 

We claim that 

w{0 > nR/io,^e{nR/io,i - 2) . (2.9.6) 
This is just the second quantized version of the inequality 

eiR/2 - d{^,xj))eiR/2 - di^,xi^^{xj))) 

> e{R/10 - di^, Xj)) ^1 - n ^(^(^' - ^/lO) j ■ 

To prove (2.9.7), we have to show that whenever xj and some Xk, k ^ 
j, are in a ball of radius R/10 centered at ^, then ball 
of radius R/2 (with the same center). Assume first that k G Jj. Then 
d{xj,x^-^) < d{xj,Xk) < R/5, whence d{^,x'^-^) < 3R/W. If, on the other 
hand. A; ^ Jj, then there exists an I G Jj such that d{xi,Xk) < R/5. Hence 
d{xj,x'^-f^{xj)) < d{xj,xi) < 2R/5, and therefore d{^,x'^^) < R/2. This 
proves (2.9.7). 

Hence, in particular, we have that 

w{0 > w{0 = w{0 e{2 - n^R/2,{) + nR/io,^e{nR/io^^ - 2)e{nsR/2,^ - 3) . 

(2.9.8) 

We now claim that 

w{^) e{2 - n^R/2,^) = nRj2,s_ {nR/2,i - l) 0{2 - n2,R/2,i) ■ (2.9.9) 
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This implies, in particular, that the operator w{$,) depends only on the Fock 
space restricted to a ball of radius 3R/2 centered at ^. Eq. (2.9.9) follows 
from the fact that if two particles with coordinates Xi and Xj are within a 
ball of radius R/2, and no other particle is in the bigger ball of radius 3R/2, 
then the two particles must be nearest neighbors. Moreover, j G Jj and 
i G Jj by construction. 

Note that (2.9.9) is a bounded operator, bounded by 2. Moreover, since 
nR/io,i ^ ^3R/2,$, we also see that 

MO-nR/io,^\<2. (2.9.10) 

Using (2.9.5), (2.9.8) and (2.9.10), we can estimate 

-2-^ / dC\\r' , -J^fev , 111- (2.9.11) 

Here wc have also used that w{^) acts non-trivially only on the Fock space 
over x-iR/2,(^- Note that the integral over the second term on the right side 
of (2.9.11) is equal to 

^deTV^[n«/io,e(n-nf)] = f IV^[N(Tf-nf)] . (2.9.12) 

Moreover, for the last term in (2.9.11), we can use (2.8.7) to estimate 

//e l|Tj,3,/,, - ^^r,x3./.Ji < 4 {hM\S{Tl^t)f" (2.9.13) 

as long as 3i? < h. 

We proceed with a lower bound on TrjF[u;(^)r2^ . In fact, we will derive 
two different lower bounds on this expression. First, neglecting the last term 
in (2.9.8) and using (2.9.9), 



(2.9.14) 
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Since is a combination of a coherent and quasi-free state, the last expres- 
sion in (2.9.14) is, in fact, easy to estimate. Let denote the one-particle 
state = E|p|<p,^pb)- We have 

= (^Tr;r[n3K/2,^Of]) + 2tr {X3R/2,i^bf + Q{'^z\iX3R/2,i^^bX3R/2,0^\^z) 
+ 3 (Tlj^[n3Ji/2,^^ll\^ (2($z|X3ii/2,^'^6X3fl/2,^|^5-) +tr(x3ij/2,5W6)^) 

<6(TT:r[n3R/2,i^(\Y • (2.9.15) 



(Here, we use the symbol tr to denote the trace over the one-particle space 
L^(A), while Tr is reserved for the trace over the Fock space.) 
A different lower bound can be obtained using 

TV^ [wiOnf] > Trjr [nR/io,e^(n^/io,e - ■ (2-9-16) 

Eq. (2.9.16) follows easily from (2.9.8) and (2.9.9). The latter trace is non- 
trivial only over the Fock space over Xr/io,^^- Denoting by IIq^ the vacuum 
on J^, we claim that 

which implies that 

^Lnno. ^ e-'^W^°)^-/3 {U{z)^U^U{z))^^^^^^^ . (2.9.18) 

Eq. (2.9.17) follows from the fact that ^b,XR/ioi ^ particle-number con- 
serving quasi-free state, whose vacuum part is given by 

exp ( - tr ln(l -|- XR/io,^^bXR/m,i)) > exp ( - tTXR/io,iUJbXR/io,i)) 

= exp ( - 47r(i?/10)^^a;/3) . (2.9.19) 

Eq. (2.9.18) implies, in particular, that 

(2.9.16) > e-^-Wio)'^-/3lV^[nK/io,^^(nK/io,^ - 2)C/(z)tn^C/(£)] . 

(2.9.20) 

The state U{z)^U^U{z) is a coherent state on Its restriction to the Fock 
space over Xr/io,(,'H is again a coherent state. In every sector of particle 
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number n, it is given by the projection onto the n-fold tensor product of the 
wave function Xr/io,^^zj appropriately normahzed. Therefore 

TTr[nR/w,Snn/w,i - 2)U{z)^U^U{z)] 

= e-{'!>z\xRno,(\'S>^) ^ (NXR/io.d^^-)" .2 9 21) 

^ n! ^ ' 

n>2 

\ z\XR/10,e z/^ J - l + {^z\XR/10.i\'^z) 

The last inequality follows from the elementary estimate x{l—e~^) > a;^/(l+ 
x) for X > 0. 

Summarizing the results of this subsection, we have shown that, for any 
< A < 1, 

Tr^[WiTf] 



+ (1 _ A)lM4e-^'^W^°)^^"/3 / ftT'^^t. ■ (2-9-22) 



The choice of A will depend on the function |$^|. If it is approximately a 
constant (in a sense to be made precise in Subsection 2.11), we will take 
A = 1, otherwise we choose A = 0. 



2.10 Interaction Energy, Part 2 

Next we are going to give an upper bound on the expectation value of W2, 
defined in (2.9.3). To start, we claim that there exists a smooth function g 
of rapid decay (faster than any polynomial) such that 

wr{x -y)< -^gidix, y)/s) . (2.10.1) 

Although WR depends on the box size L, g can be chosen independent of 
L for large L. This follows immediately from the following considerations. 
First of all, we have, from the definition (2.6.2) of Jr and because of i? < s, 

fR{x)<R sup |V%)|. (2.10.2) 

d{x,y)<s 
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Recall that h{x) = \A\-^ Ep(l-^^(«P))e~'^^, where 1 — is a smooth function 
supported in a ball of radius 2. We need the following elementary lemma. 

Lemma 5. Let o : — C 6e a smooth function, supported in a cube of side 
length 4, and let u{x) = \A\~^Y^p0{sp)e'~'^P^. Then, for any non-negative 
integer n. 

Here, A denotes the Laplacian on M^, not on A. 
Proof. Introducing coordinates x = (.xi, .X2, X3), we can write 
/ 3 \ " 



—ipx 



u{x) 2L^ ^(1 - cos(27rx,/L)) = \K\-^ ^(-Ad)"o(sp)e 

\ i=l ) V 

(2.10.4) 

Here, A(i denotes the discrete Laplacian in momentum space, which acts as 
L^'^{-Ad}f{p) = 8f{p) - E|e|=i f(P + 27re/L). It is easy to see that the 
function (—Ad)"/ is bounded by ||(— A)"/||oo- Moreover, if / has support 
in a cube of side length i, then (— A(j)"/ is supported in a cube of length 
£ + 47rra/L. This implies that 

1(2.10.4)1 <s2"||(-A)«o|U f- + 2^V . (2.10.5) 

\Trs Li J 

On the other hand, note that 1 — cos{2TTXi/L) > 8L~^ m.mk^z\xi — kL\^, 
and hence 

3 

2L2 ^(1 - cos(27ra;i/L)) > 16 d{x, Of . (2.10.6) 

i=l 

This proves the lemma. ■ 

Applying the lemma to the function V/i in (2.10.2), and using the defi- 
nition (2.6.3) ofwR, we immediately conclude (2.10.1). 

We now decompose the function g into an integral over characteristic 
functions of balls. Such decompositions have been studied in detail in [3]. 
Recall that j is defined in (2.7.1). According to [3, Thm. 1], we can write 



9{t) = / 
Jo 



drm{r)j{t/r), (2.10.7) 
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where 

m{r) = ^r{g"{r)-rg"'{r)) . (2.10.8) 

Note that m is a smooth function of rapid decay. Since j is monotone 
decreasing, we can estimate 

git) < jit) J dr\m{r)\+J dr \mir)\ jit/r) . (2.10.9) 

This estimate, together with (2.10.1), impUes that 

W2 < — ^ r dr (5{r - s)J^ dt |m(t)| + V(^A)l) 

„ oo n 

^ / 0E E Xr/2d^j)Xr/2,^iXi) ■ (2.10.10) 
•^^ n=0 j=l ieJj 

Let VriO denote the integrand in the last Unc in (2.10.10). Because 
dixi,Xk) > R/5 ior i,k £ Jj, the number of ,Xj inside a bah of radius r/2 is 
bounded from above by (1 + 5r/R)^. Thus we have 

VriO<nr/2A^ + 5r/R)\ (2.10.11) 

Moreover, we triviahy have that Vri^) < Wr/2,$('^r/2,{ ~ !)• By combining 
these two bounds, we obtain 

VriO < /K/2,«) , (2.10.12) 

where 

f{n) = nmin{n - 1, (1 + 5r/R)^} . (2.10.13) 

Proceeding similarly to (2.9.11), using that |/(ra) — ri(l + 5r/R)^\ < (1 + 
(1 + 5r/i?)^)/4, we can estimate 

1V^K(e)n] <Tr^[/(n./2,c)Tg 

< Trr[finr/2,i)nt] + (1 + 5r/R)^Tr^[nr/2A^i - f^f)] 
+ i (1 + (1 + 5r/R)')' ||Tj^^^^^^ - nl^^^Jl, . (2.10.14) 

When integrating over ^, the last two terms can be handled in the same way 
as in the previous subsection, see Eqs. (2.9.12) and (2.9.13). We have to 
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assume that r <b, however. For the first term on the right side of (2.10.14), 
we estimate 

< min {TTr[nr/2,iinr/2,i - 1)^11] , (1 + 5r/i?)3Tr^[n,/2,^^^^l } . 

(2.10.15) 

Similarly to (2.9.15), 

Tr^[nr/2,i{nr/2,i - mt] < 2 (Tr^[n,/2,fO^^)' , (2.10.16) 
and hence 

Moreover, 

Tr^[nr/2,^nt\ = ^r^e^ + {^^\Xr/2,^\^z) . (2.10.18) 
Using convexity of the function x x'^ /{I + x), we obtain the bound 

(2.10.19) 

We use (2.10.19) in (2.10.14) and integrate over ^. We obtain (assuming 
r < 6, as mentioned above) 

JC'^rivrmi] < {l + 5r/Rf^r'TT^[N{Ti-ni)]+8\A\ (Ir^Qu^f 



+ 



(l + (l + 5r/i?)3)'(63|A|5(T^,Of)) 



1/2 



In order to be able to compare the last term with the last term in (2.9.22), 
we note that 

< f 1 + 51 / daxR/io,^+a ■ (2.10.21) 

V -n-/ J\a\<r/2+R/10 



30 



Here, we denote by f the normalized integral, i.e., we divide by the volume 
of the ball of radius r /2 + R/W. Using monotonicity and convexity of the 
map X I— x"^/ {1 + x), we thus have the upper bound 



l + 4{^^\Xr/2,^\^z)/il + 5r/Rr 

<ri + ^V/ ^^J^dx^A0,i±^^_ ^2.10.22) 

V RJ J\a\<r/2+R/10 ^ + M^z\XR/10,^+a\^z) 

After integration over ^, the right side of (2.10.22) simply becomes 

5^\6 . {<^s\xr/w,^\'^z)^ / /^6r\^ /• {'^z\XR/io,i;\^z)^ 



Rj Ja ^ l + A{^g\xR/lo,i\^z) ~ \ RJ J a 'i- + {^z\XR/10,i\^z)' 

(2.10.23) 

Here we have used the fact that r > s > R for the relevant values of r. 

As noted above, the estimates leading to (2.10.20) arc only valid for 
r < b. To bound the expectation value of W2 in (2.10.10) we have to 
consider all r > s, however. For r > 6, we use (2.10.11) to obtain the simple 
estimate 

The contribution of r > 6 to the integral in (2.10.10) is thus bounded from 
above by 



- r dr\m{r/s)\ [ deiV^K(6n] 
S Jb J A 



- I'' drr''\m{r)\ . (2.10.25) 

Since \m\ is a function that decays faster than any polynomial, the last 
integral is bounded above by any power of the (small) parameter s/h. 
Let c denote the constant 

pi POO 

c = J dr\m{r)\ + J drr^\m{r)\. (2.10.26) 
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+ 



To summarize, we have derived the upper bound 

+ 6^— ^TV^[NJ^f] / drr^\m{r)\ 

S^i^S Jb/s 

+ 8-^^ / . (2.10.27) 

2.11 Interaction Energy, Pcirt 3 

We now put the bounds of the previous two subsections together in order to 
obtain our final lower bound on Trjr[WT^]. We will distinguish two cases, 
depending on the value of a certain function of |$^|, given in (2.11.1) below. 
Assume first that 



' ^.^v.v..,,, ^ _|A|(ii;'^^)^ . (2.11.1) 

This condition means, essentially, that |$^| is far from being a constant. In 
this case, we choose A = in (2.9.22), and find that the contribution of the 
last terms in (2.9.22) and (2.10.27), respectively, is bounded from below by 

8.|A|a,^(^^--(-j ..-8c^j. (2.11.2) 

Next, consider the case when (2.11.1) is false. In this case, using (2.10.21) 
for r = 3i?, as well as convexity of x i-^ a;^/(l + x), we find that 

Pick some > 0, and let H C A denote the set 

H = G A : {^,\x-m/2,i\^z) > IQ^DB?q} . (2.11.4) 
Using (2.11.3), as well as monotonicity of x i— > a;/(l + x), we find that 

JJ^{^z\X3R/2,i\^z) < —-\A\R'Q{l + DR'e) . (2.11.5) 
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Similarly, we have the estimate 

\B\<\K\^'^{l + DB?e) . (2.11.6) 
We choose A = 1 in (2.9.22) and estimate the relevant term from below by 

> j^^^ (Tr^K2,c(nR/2,^ - ml] " 6 (Tr^[n3^/2,4l^f])') . (2.11.7) 

Using Tr^[n3K/2,5l7f] = 97ri?3^^/2 + {<^^\x3R/2,i\'^z), the definition of B in 
(2.11.4) and convexity of x i— > x^, we can bound the last term as 

J^^d^ (TT^[nsR/2,^nl]y < 4|A| (^Y^'^^) + 187r|f|2/?3 ^iQ^DR^gf . 

(2.11.8) 

We now investigate the first term on the right side of (2.11.7). A simple 
calculation shows that 

T^r[nR/2,d'^R/2,^ - l)f^f] 

= T^T[nR/2,^{nR/2^(. - l)nh] + 2{<^^\xR/2,^u)bXR/2,^\^z) 

+ ^R^Q. {'^z\XR/2,i\^z) + {'fz\XR/2,^\'^zf ■ (2.11.9) 

Here, we have used again the translation invariance of Qh- Note that this 
invariance also implies that the first term on the right side of (2.11.9) is 
independent of ^. Since is a quasi free state, it can be rewritten in terms 
of the one-particle density matrix ujb as 

T^T[nR/2,^{nR/2,^-l)nb] = [tr XR/2,^i^bf +tT XR/2,^i^bXR/2,^^b ■ (2.11.10) 

The first term is just {irR^g^/G)'^, and the second is bounded from above by 
this expression. Therefore, 

j^diTTr[nR,2AnR,2,i - l)^b] < m {iR'^Qo.y ■ (2.11.11) 

Note also that g\xR/2,^^^bXR/2,^\^ g) < tTXR/2,^^^b {^z\Xr/2,^\^z), and thus 

['^{'^z\XR/2,i'^bXR/2,i\'^z) + ^R^Q^ {^z\XR/2,i\^z)) 
< "^R^Qu. di {^s\XR/2,e.\^z) ■ (2.11.12) 
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The last expression can be bounded from above using (2.11.5). For the 
last term in (2.11.9), we use Schwarz's inequality, together with (2.11.5), to 
estimate 

J^^^ {^z\XR/2A^zf > ^ {^z\XR/2,i\^z)^ (2.11.13) 

Here we set again = \z\'^/\A\. 

Putting all these estimates together, we have thus derived the lower 
bound 



Ja\b 

- |A|— ^^c^^Y^ (1 + DR^q) +^j^dC {^^\XR/2,i'^bXR/2 



The first integral on the right side of (2.11.14) can be rewritten as 



|A|- 



144 

where we introduced the notation 



jjx\ub{x)\'j{dix,0)/R) > |A|^76', (2.11.15) 



76 = ^ Jjxu;b{x)j{d{x,0)/R) . (2.11.16) 

Eq. (2.11.15) follows by applying Schwarz's inequality to the integration over 
A, noting that Jj^dx j{d{x,0)/R) = AttR^. 

It remains to integrate the last term in (2.11.14). We claim that 

f tt'^R^ 
d^ {^z\XR/2,i^bXR/2,^\^z) > 1^'^^ (76 - RPcQu^) ■ (2.11.17) 
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To see this, we write 

^ J^(^{^z\XR/2,^^^bXR/2,^\^z) - \z\^ j^dxujb{x)j{d{x,{))/R) 

= I dxdy {^g{x + y)* - ^M*) ^M^b{x)j{d{x, 0)/i?) 
^AxA 

> -W^sh [ dx \\^z{x + ■)- ^z{-)\\2H{x)\j{d{x,Q)/R) . (2.11.18) 
J A 

We can estimate |wb(x)| < a;(,(0) = q^j. Moreover, writing the 2-norm as 
a sum in momentum space, and using the fact that $5- has non-vanishing 
Fourier coefficients only for |p| < pc-, it is easy to see that ||^*^(x + •) — 
*J'z(")ll2 ^ ||*J'z||2PcC^(a;, 0). Since the range of j( • /R) is R, the integral over 
A can be estimated as jj^dx j{d{x,0)/R)d{x,0) < R J^dxj{d{x,0)/R) = 
AttR"^. This yields (2.11.17). 

Collecting all the estimates above, we conclude the following lower bound 
on the expectation value of W: 



TV^[WT;] > 24-^IV:r[N(Ti - f!f)] - 6h 

-Hi|5(63|A|S(n.nf))"^(8 + (l+63fc|; 

-«|A|(8£c|:4^£...V«l) 

+ 47ra'|A| min {.Ai , A2} ■ (2.11.19) 
Here we have used the simple bound a' < a, and we have set 
nf Att / R\^ a6^R^\ 

= (10) ""-'^VlJ^j <'-"-'°> 

and 

■A2 = {qI + 2gsjb + 7b) + '^QujQzi'i- - RPc) 

+ ^2 ^1 - -^y (1 + DR^g) - 2Tr 3^ R^Q^^ 

-Qze i-^{l6^DyR^Q+——{l + DR^Q)j . (2.11.21) 



/la' . R 
(\^r.— 
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We will choose D = {Epg)~^/^ in order to minimize the error terms in A2- 
Moreover, since a' /a contains a factor (1— e) (sec (2.7.7)), it is best to choose 
e = R/s. Wc note that one can also use the simple bound (2.4.10) in order 
to estimate in the error terms. 

Since Rq <^ R <^ s, the term in round brackets in the first line of 
(2.11.19) is non-negative and, therefore, we will need a lower bound on 
TrjF[N(T| — ^f)] ■ Moreover, we will need an upper bound on the relative 
entropy >S'(T^,0^). Appropriate bounds will be derived in the next two 
subsections. 

2.12 A Bound on the Number of Pcirticles 

Our lower bound on the expectation value of W in the previous subsection 
contains the expression TrjF[N(T^ — ^^f)] , multiplied by a positive parame- 
ter. Hence we need a lower bound on this expression in order to complete our 
bound. In fact, we will combine the first term on the right side of (2.11.19) 
with the last term |Trjr[ICT^] in (2.7.15), which we have not used so far. 
I.e., we seek a lower bound on 

24^Tr^[N(n - m - ^'4) + ^Tr^[ (N " N)^ . 

(2.12.1) 

(Here we have used that e = R/s, as mentioned at the end of the previous 
subsection.) First, note that TVjr [NQf] = |z|2-FTrjr[Nfi6] = |zp-FTr;r[N07r] 
and Tr^[NT|] = |zp + Tr^[NT^]. Let N> = E|p|>p,4ap denote the 
number operator on J^^. Using that fi^ = H (g) r° and that = H (g) F^, 
we can thus write 

'iv^[N(Tf - nl)] = Tv^^ [N>(r,-- r°)] . (2.12.2) 

For the second term in (2.12.1), we use 

(N - > {\z\^ + Tr^^ [N>r°] - iV)^ (2.12.3) 
+ 2 {\z\'^ + TV^> [N>r°] - iV) (N - \z\^ - TV^> [N>r°] ) , 

and hence 

Tr^[(N - NfT^ > + Tr^> [N>rO] - N)^ (2.12.4) 
+ 2 {\zf + Tr^> [N>rO] - N) Tr^^ [N> (F,- - F^) ] . 
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Thus, we conclude that the expression (2.12.1) is bounded from below by 
27raC ^1 _,|2 ^ ^^>j.o] _ j^f + t^^^ [(n> _ ^0)p^^] (2.12.5) 



|A| _ 

24 / a' o aR\ AnaC ^ r^.^^m 



X 



We will choose i?, s and C below in such a way that R^g <^ 1/C and 

i? <C s. The last term in square brackets is thus positive, irrespective of the 
value of 1^1 . Hence we need to derive a lower bound on Trjr> [N-^ (Fg — F^) ] . 
To this end, we note that, for any /li < 0, 

S{T,, rO) - /?/xTV^> [N>r,-] > /3(/(m) - 7(0)) . (2.12.6) 

Here, we denoted 

/(/^) = E (l - e-^(^'-''°-'^)) . (2.12.7) 

|p|>Pc 

Eq. (2.12.6) is simply the variational principle for the free energy /(//). Since 
/ is completely monotone (i.e., all derivatives are negative), we can estimate 
m > 7(0) + ^7(0) + ^^^"(0). But 7(0) = -TV^> [N>rO] and hence, 
optimizing over all (negative) /j,, 

iv.,r(r,-r»)] l5(r.r») E ' 

\ \P\>PC , 

(2.12.8) 

We can use (2.4.6) to estimate the relative entropy as ^(r^, F^) < 8-7r|A|a/3£)^. 

For the sum over p, we use that cosh(x) — 1 > jl. In the thermody- 
namic limit, we can replace the sum over p by an integral. We thus have to 
bound 

We use two different bounds. On the one hand. 



1 /■ , 1 vr^ 



^T?T-2 yi=^=- (2.12.10) 

On the other hand, since /xq < 0, 

f I r 1 An 

/ dp-^ -2</ dp-^ = — . (2.12.11) 

Ap\>pc ip -f^or J\p\>po ip ) Pc 
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In combination, these bounds imply that 

(2.12.9) < ^ ( . / . ^ . (2.12.12) 

Using (2.12.12) in (2.12.8), we obtain the lower bound 

IV^> [N> (r,— rO) ] > -const. |A|^ {p, + y^)"'/' - o(|A|) . 

(2.12.13) 

We apply this bound in (2.12.5), noting again that the last term in square 
brackets is positive. We conclude that 

(2.12.1) > {\z\^ + TV^> [N>r°] - Nf - - o(|A|) , (2.12.14) 

with 

^3 /2 r 1 

= const. |A|^ {pc + V^y^'^ -^ + c[q(i + 2/Vc) + Qu) 

(2.12.15) 

Here we have used (2.4.10) to bound \z\^ from above, as well as the fact that 
2.13 Relative Entropy, Effect of Cutoff 

Next, we are going to give an estimate on the relative entropy of the two 

states T| and Vl^. This is needed for the lower bound on the expectation 
value of W obtained in (2.11.19). As already noted, the relative entropy is 
invariant under unitary transformations, and hence 

SO:i,Q.i)=S{Tl®Vg,^j,). (2.13.1) 

We want to bound this expression by the relative entropy of 11 (g) with 
respect to = n (g) r°, which satisfies 

5(n ® r^, n^) = S{V^, r°) < 87ra/?|A|^2 (2.13.2) 

according to (2.4.6). I.e., we want to estimate the effect of the cutoff h on the 
relative entropy 5(11 ® F^, Q;,). Here it will be important that H is not the 
vacuum state. The cutoff h corresponds to a mollifying of the one-particle 
density matrix in momentum space, and the error in doing so would not be 
small enough if this one-particle density matrix is strictly zero for \p\ < pc- 
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This is the reason for replacing the vacuum state Hq by a more general 
quasi-free state 11 in Subsection 2.5. 

If r2(^ denotes a general (particle number conserving) quasi-free state 
with one-particle density matrix w, it is easy to see that S'(r,Otj) is convex 
in Lo for an arbitrary state F. The one-particle density matrix of can be 
written as 

^'^ = ■r^Xl'?*(^)lZ 2('^'r(P + g) +'^7r(p- g))|p)(p| • (2.13.3) 
' ' q p 

Therefore, 

s{u ® r,^, f],) < i- ^ r}b(.q)S(n r,-., n^) , (2.13.4) 

where Qq is the quasi-free state corresponding to the one-particle density 
matrix with eigenvalues ^{u>Tj-{p + q) +u!Tr{p — Q))- (This is the same estimate 
as in [17, Sect. 5.2].) Moreover, as in [17, Eq. (5.14)], simple convexity 
arguments yield 

s{u rg,nq) < (1 + 1-^) s{u r^, n^) (2.13.5) 
p ^ 

for any t > 0. Here 

= 1, 2 + u {p + q)+u {p- q) 

Note that ho{p) = £{p), defined in (2.8.1). To estimate the expression 
(2.13.5) from above, we need the following lemma. 

Lemma 6. Let £ : M.^ M+, and let L± = ± sup^ sup^ n^n^i ib(Q'V)^£(p) 

denote the supremum (infimum) of the largest (lowest) eigenvalue of the 
Hessian of i. Let lo-jt{p) = [e^^P^ — 1]~^, and let hq{p) he given as in (2.13.6). 
Then 

hq{p) - ho{p) < L+q^ , (2.13.7) 



and 



hqip) - ho{p) > q^L- + q' min{L_, 0} - 4g^ sup[\Viip)\^co^ip)] 

p 

- 2q^{\q\ + |p|)2sup[|V£(p)|Vp'] • (2-13.8) 
p 
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Proof. By convexity of x i— ln(l + 1/x), 

h,{p) < '^{e{p + q)+ e{p - q)) < i{p) + L+q^ , 



(2.13.9) 



proving the (2.13.7). To obtain the lower bound in (2.13.8), we proceed 
similarly to [17, Lemma 5.2]. We can write 

hg{p) - hoip) = ^ dX{l- X)^hxgip) . (2.13.10) 

Denoting p± = pzLXq and uj± = a;,r(p±)) we can write the second derivative 
of hxq{p) as 



■^hxqip) 



1 



d 



_(w++a;_)2 (2+a;+ + a;_)^J V^^ 

2 

(lj+ + a;_)(2 +LJ++ a;_) dX^ 



^(c.+ + ^- 
(w+ + a;_). (2.13.11) 



The first term is positive and can thus be neglected for a lower bound. For 
the second term, we use 



^a;+ =u+{l+ u+) [(1 + 2u+) {qVi+f - {qV)H, 



(2.13.12) 



where we denoted ^_|_ = (.{p + \q). The last term in the square brackets is 
bounded above by —q^L^. Moreover, 

(1+2^+) {qVi+f 

<g'(|p| + M)'sup[|V£(p)|Vp'] +2g'sup[u;,(p)|V£(p)|2] . (2.13.13) 
p p 

The same bounds hold with + replace by — . Using in addition that 
1 a;+(l + a;+) + a;_(l + (J 



- < 



2- (a;+ +a;_)(2 + a;+ +a;_) 



< 1, 



(2.13.14) 



we arrive at (2.13.8). 



Let g : t-^ [0, 1] be a smooth radial function supported in the ball of 
radius 1. We assume that g{p) > \ for \p\ < |. We then choose 



Hp) = PiP^ - ^o) + Pplgip/Pc) . 



(2.13.15) 



Since i{p) = ln(l + 1/%) for \p\ < pc by definition, this corresponds to the 
choice TTp = (exp(/3(p2 — ^q) + PPcdip/Pc)) ~ 1)^^- Note that, in particular. 
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TTp < const. and hence P < const. M/{Ppl) ~ Pc\A.\/(3. This bound is 

important for estimating the error term Z^^^ in (2.5.6). 

For the £ given in (2.13.15), both p-^L+ and p'^L^ are bounded in- 
dependent of all parameters. Moreover |V^(p)| < const. Using that 
^7r{p) < ^{p)~^ < {Pp^)~^7 the bounds in Lemma 6 imply that 

-DPq^ (1 + f5{\p\ + \q\f) < hgip) - ho{p) < D(5q^ (2.13.16) 

for some constant D > 0. 

Using the fact that sinh(a;)/a; < cosh(x) for all x G M, we can estimate 

{h,(jp) - hoip)) - ^HM^Ti) (2.13.17) 

The estimate (2.13.16) implies the bound 

{hM - hoip)f < D\Pq^f (1 + P{\p\ + \q\ff . (2.13.18) 

Moreover, using the upper bound in (2.13.16) to estimate hq{p) — /to(p) from 
above, we obtain 

= u\p) + i {oj^ip + q)+ u^{p - q)) (1 + 2j{p)) (2.13.19) 

as an upper bound to the last fraction in (2.13.17). Here, we denoted a;*(p) = 
\^Qho{p)-D(}tq assuming t to be small enough such that ho{p)—Dptq'^ > 

for all p. Recall that /io(p) = ^ip) is given in (2.13.15). 

In the thermodynamic limit, the sum over p in (2.13.5) converges to the 
corresponding integral and, hence, we are left with bounding 

/ dp {i+m + wiff {^\p) + UMp+9)+Mp-<i)) (1+2c^*(p))) 

(2.13.20) 

from above. We can replace ci;,r(p — q) by a'7r(p + q) in (2.13.20) without 
changing the value of the integral. Using Schwarz's inequality, the fact that 
(j^nip) < i^^{p), and changing variables p ^ p — q, we see that (2.13.20) is 
bounded from above by 

2/ dp{l + P{\p\+2\q\f)^Jip){l + 2Jip)). (2.13.21) 
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It remains to bound a;*(p). For this purpose, we need a bound on i{p) 
Dptq^ for an appropriate choice of the parameter t. We will choose t ■ 
min{l, {b^q^)~^}- We then have tq^ < and it is easy to see that 



iip) - D(5tq^ > 13 



(2.13.22) 



in this case. Since l{p) > P{p'^ — no), this can be seen immediately in the 
case \p\ > ipc- For \p\ < pc/2 even a slightly better bound holds, this time 
using the fact that g{p) > 1/2 for \p\ < 1/2. 

We will choose b and Pc below in such that a way that bpc S> 1. Denoting 
by r the (positive) number 

r = -/3/xo + /3p^Q-^) , (2.13.23) 



we thus have the bound 



u'ip) < 



(Je\^P^ - 1 



< ( 1 + \ 1 . (2.13.24) 



2f 

The last bound follows from the elementary inequality (e^ — 1)~^ < e~^(l + 
1/x) for all X > 0. We insert this bound for a;* into (2.13.21). Simple 
estimates then yield 

(2.13.21) < const. ^ (} + r'^^^) (l + iPl^?) ■ (2.13.25) 

Combining (2.13.5), (2.13.17), (2.13.18) and (2.13.25), and using that 
we have thus shown that 

5(n ® Tg, Qg) < (2 + h^q^) s{n ® Ts, 

+ const. |A|/3^/\V-^/2 (1 + {Pq^f) + o(|A|). (2.13.26) 

After inserting this estimate in (2.13.4) and summing over q, this yields the 
bound 

5(Tf , Qf) < const. |A| [aj^Q^ + ^^^^-^^ + o(|A|). (2.13.27) 

Here, we have used the assumed smoothness of rj to estimate %('7)|'7|" < 
const. |A|6~" for integers < 8. We have also used the fact that we will 
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choose b"^ f3 and hence, in particular, < const. Moreover, the as- 

sumption (2.4.6) has been used to bound S{Il Tg, = S{Tg, T^). 

We have thus shown that the effect of the cutoff b on the relative entropy 
can be estimated by |A|(/3/r)^/^5~'^. We note that the exponent —4 of b 
is important, since the relative entropy has to be multiplied by b^ in the 
estimate (2.11.19). 



2.14 Final Lower Bound on Fs{f3) 

We have now obtained all the necessary estimates to complete our lower 
bound on Fg{(3). For this purpose, we put all the bounds from Subsec- 
tions 2.7, 2.11, 2.12 and 2.13 together. In fact, from (2.7.15), (2.11.19), 
(2.12.14) and (2.13.27), we have the following lower bound on Fs{f3): 

Fm > - InTV^^ exp ( - /3T^(z)) - Z^^^ - Z^^^ - Z^^) 
-(x-xO p%-o(|A|) 

|P|<J'C 

+ ^(kT + iv^>[N>rV^)' 

+ 47rS|A| min {2^^ ^ ^2 ^ 2^.(^^ + q^) + ^2 ^ ^ ^2.14.1) 
where we denoted 



Z^'') = const. a|A| 



1/2 



(2.14.2) 



Here, we have used the definition (2.7.7) of a', (2.4.10) to bound Qg in the 
error terms, together with 7^ < and lim|A|_»oo Quj ^ Q- This last estimate 
follows from (.{p) > /3(p^ — ;Uo) and (1.6). The error terms Z(^) and Z(^) are 
defined in (2.5.6) and (2.12.15), respectively. 

Using (2.7.8), the term in the second line of (2.14.1) can be estimated 
by (x- x')E|p|<p,p'7rp < (x- x')p2p < const. {a/R)'^pl\A\/^. Here, we 
have also used the bound on P = % derived after Eq. (2.13.15). 

The two terms in the third and forth line of (2.14.1) can be bounded 
from below independently of z, simply using Schwarz's inequality. More 
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precisely, introducing = |A| ^Trjr^ [N>r°] = Qu — -P/|A|, we obtain 

+ Tr^, [N>rO] -N)' + 4^a\A\ {gl + 2g,{j, + g^) + gl + j^) 



|A| ^ 

(2.14.3) 



We note that 

^0 



^" = 7^3/ ,fa._L) 1 + (2-14.4) 

71^)"^ J\p\>Pc et='(P t'") - 1 



in the thermodynamic hmit. Hence, from (1.6), 

= min {, , ,,(T)} - ^ 1^^^^^ dp ^^pz^ + o{l) 

> minie , g,iT)} - + o(l) . (2.14.5) 

This estimate is obtained by using e^^'^~^°^ — 1 > l^p^ in the denominator 
of the integrand. Moreover, g^ < g^, < min{^, gciT)} + o(l). 

It remains to give a lower bound on 7;,. According to (2.11.16) and 
(2.8.2), 

lb = ^ dx u^ixMdix, 0)/b)jidix, 0)/R) . (2.14.6) 
We note that, since UTr{x) is real, 

(2.14.7) 

We can estimate d{x,0) < Rm the integrand in (2.14.6). Since £{p) > 

Ir/I < 1 and Jj^dxj{d{x,0)/R) = AttR^, the contribution of the last term to 

(2.14.6) is bounded by 

in the thermodynamic limit. Moreover, we can bound 77 from below as 
vii) > 1 ~ const, i^, and hence 



Quj>lb> 



( R^\ 
g^jil- const, -p- j - const. R^p'^^'^ - o(l) . (2.14.9) 
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Using the bounds on and 76 just derived, we have 

(2.14.3) > 47ra|A| (2q^ -[q- Qc{T)]\) (2.14.10) 



- const. S|A| [g^ + ^) + (| + ) " • 

In particular, the terms in the third and forth line of (2.14.1) are bounded 
from below by the right side of (2.14.10). 

2.15 The "Free" Free Energy 

We now insert the lower bound on F^{f3) derived in the previous subsection 
into (2.3.9). We note that the only ^-dependence left is in the first term in 
(2.14.1), which is the "free" part of the free energy. Taking also the constant 
PLqN in (2.3.9) into account, we are thus left with evaluating 

- ^ In / d'^z TV^> exp ( - /3T,^(z)) 
= MoiV + i I 1; In (/3e(p)) + ^ In (l - e'^'^^^) j . (2.15.1) 

\|p|<Pc |p|>Pc / 

Using X > (1 — e~^) for non-negative x, (2.15.1) becomes, in the thermody- 
namic limit, 

(2.15.1) > iV^o + / dp In (1 - e-^^(f)) - o(|A|) . (2.15.2) 

Recall that e{p) is defined in (2.7.10). It is given by e{p) = (1 — x+x')p^— /xo 
for \p\ < 1/s, and satisfies the bound e{p) > x'p^ for \p\ > 1/s. Hence 

/ dp In f 1 - e-l^^pA > (1 - X + x')~^/^ / dp In f 1 - e-^(p'-^°)) 

(2.15.3) 

We will choose <C x'/? below, in which case the last integral is expo- 
nentially small in the parameter s^/(x'/3). Inserting the definition (1.4) of 

/o(/3, q), we have thus shown that 

(2.15.1) > |A|(1 - x + ^'y'^/^foif3, q) (2.15.4) 
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Here, we have also used that fio < 0. 
2.16 Choice of Pcirameters 

We have now essentially finished our lower bound on g). It remains to 
collect all the error terms, and choose the various parameters in an appro- 
priate way. All the error terms we have to take into account are given in 
(2.1.1), (2.3.9), (2.14.1), (2.14.10) and (2.15.4). 

We will choose the various parameters in our estimates as follows: 

i?=£.-V3(„^2^5/2)3/403 6 = ;3l/2(„,2^5/2) -121/403^ 

S = (/3^-V3) 1/3 (^^2^5/2) 1/403 _ ^^.16.1) 

Moreover, 

ip = a{ag^P'/y^ , C = {ag^p'/y (2.16.2) 

for 4/403 < A< 79/403 and 2/403 < B < 161/403. Depending on fiQ, we 
choose 



Pc 



rV^(a,^/?^/^)«V403 ^i^^i < (,^2^5/2)162/403 ^^ ^^ ^^ 

otherwise . 



Finally, we choose x = ^ {ag'^j3^/'^) ^ for some 5 > Q. Our estimates 
then imply that 

/(/5, q) > /o(/9, q) + 47ra {2g' -[g- e,{(3)]l) (l - o(l)) , (2.16.4) 

with 

0(1) < Q((/?£.2/3)-l)(a^2^5/2)2/403-. ^2.16.5) 

for some function Cs, depending on 5, that is uniformly bounded on bounded 
intervals. 

The choice of the parameters Pc, b, s, R and k is determined by mini- 
mizing the sum of all the error terms. The main terms to consider arc, in 
fact, the terms Mp^ ~ \A\pl in Z^^^ in (2.3.7), and \A\ag'^{:K + R/ s) as well 
as \A\aR-^{b^al3g'^ + (pi - 110)-^/%-^)^/'^ in Z^^^ in (2.14.2). Moreover, we 
have to take the restriction <C x/3 in (2.15.3) into account. This leads to 
the choice of parameters above. 

All other error terms are of lower order for small afp'fi^l'^ . This is true, in 
particular, for all the terms containing and C, which explains the freedom 
in their choice above. 
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2.17 Uniformity in the Temperature 

Our final result, Eq. (2.16.4), does not have the desired uniformity in the 
temperature. It is only useful in case the dimensionless parameter ag^P^I'^ is 
small. In particular, one can not take the zero-temperature limit Pg'^^^ oo. 
The reason for this restriction is that our argument was essentially per- 
turbative, using that the correction term we want to prove is small com- 
pared to the main term, i.e., ag'^ ^ fo{P,Q)- Below the critical temper- 
ature, fo{l3,g) = const. hence the assumption is only satisfied if 

If the temperature is smaller, we can use a different argument for a lower 
bound on /, which uses in an essential way the result in [11]. There, a lower 
bound in the zero temperature case was derived. 

To obtain the desired bound for very low temperature, it is possible 
to skip steps 1-5 entirely, and start immediately with the Dyson lemma. 
Lemma 2, applied to the original potential v. Using this lemma, we have 
that 



Hn>Y1 



N 

Vjil-il-H)x{Pjf)Vj (2.17.1) 



+ {1-£){1- x)aUR{d{xj,x^^{xj))) - - ^ WR{xj - Xi) 



Since d{xi,Xk) > R/5 for i,k e Jj, we can estimate (using (2.10.1)) 

TV 

Ea \ ^ aN 
- 2^ WR{xj - Xi) < const, . (2.17.2) 

j=l i&Jj 

Moreover, the calculation in [11] shows that, for the choice x = (a^£»)^/^^ 
and R = a{a?g)-^/^\ 

N 

^ --Aj + (1 - e){l - ><)aUR{d{xj,x-^^{xj))) 

> 47raA^e(^l -£- const, (a^e)^/^^) . (2.17.3) 

(Strictly speaking, this result was derived in [11] for Neumann boundary 
conditions, and with Jj = {1,...,A'^} independent of all the particle co- 
ordinates. It is easy to see, however, that the same result applies to our 
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Hamiltonian, being defined with periodic boundary conditions, and having 

a shghtly smaller interaction.) 

For this choice of x and R, we thus have 

i?iv > ^ li^/^A]) + AnaNg ( 1 - e - const, {a'^g)^^^'^ - const. J^^j > 

(2.17.4) 

with l{\p\) = — x/2 — (1 — x)x(p)^)- To obtain a lower bound on the 
free energy for this Hamiltonian, we can go to grand-canonical ensemble, in- 
troducing a chemical potential in the usual way. Taking the thermodynamic 
limit, this yields 



f{(5, g) > sup \f,g + [ dp In (l - g-^Wf)"'^)) | 

-I- Airag^ ^1 - £ - const, (a^g)^^^'^^ - const, -^{a^gf^^'^ . 

(2.17.5) 



Recall that xip) = i^i^p), where v is a function with < iy{p) < 1 that is 
supported outside the ball of radius 1. This implies that l{p) = (1 — ^>i)p'^ 
for \p\ < 1/s, and l{p) > ^xp^ for \p\ > 1/s. Hence 

/ dp In f 1 - e-^W^')-'^)) > (1 - ix)-3/2 / dp In f 1 - e-^^^"-^A 

(x/3)"^/^ J\p\^>^p/s'^ \ 

(2.17.6) 

The last expression is exponentially small in the (small) parameter s^/(x/3). 
We choose, for some (5 > 0, 

- Bia^o)^'^'^^^ - ^^^^^^'^'^ (2 17 7) 

s -p[a g) , e - ^^^2^5/2)2/5' ^^■^'■0 

and obtain 

/(/?, q) > MP, q) + ATrag'' (1 - o(l)) , (2.17.8) 

with 

»(1) = const. (1 + -^L.) ^ ^-0^) . (2.17.9) 
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Compared with our desired lower bound, we also have to take into account 
the missing term 0(0^ — [g — gc{P)]\) , which can be bounded as 

a(^^ -[g- gc{f3)]l) < const, agp-^^"^ . (2.17.10) 

In combination, the estimates (2.16.4) and (2.17.8) provide the desired 
uniform lower bound on the free energy. Depending on the value of 
one can apply either one of them. To minimize the error, one has to apply 
(2.16.4) for a^?2/55/2 < (^3^) 403/6885^ (2.17.8) otherwise. This yields our 
main result. Theorem 1, for a = 2/2295 — 6. 

A Appendix: Proof of Lemma 2 

For simplicity, we drop the ~ on v and a in our notation in this appendix. 
We start by dividing up A into Voronoi cells 

Bj = {xeA : d{x, Vj) < d{x, yk) VA; ^ j} . (A.l) 

For a given G iJ^(A), let ^ be the function with Fourier coefficients ^{p) = 
x{p)'^{p)- We thus have to show that 

/ dx[\Vax)f + lv{d{x,yj)mx)f]>{l-e)a[ dxU{d{x,yj))mx)f 
JBj JBj 

-- dxwR{x-yj)\'ip{x)f . 
£ Ja 

(A.2) 

The statement of the lemma then follows immediately by summing over j 
and using the positivity of v. 

We will actually show that (A.2) holds even when the integration region 
Bj on the left side of the inequality is replaced by the smaller set Br = 
Bj D {x e A : d{x,yj) < R}. Note that the first integral on the right side 
of (A.2) is also over this region, since the range of U is supposed to be less 
than R. 

As in Subsection 1.3, let cpy denote the solution to the zero-energy scat- 
tering equation 

-A^x) + ^v{\x\)Mx) = (A.3) 

subject to the boundary condition lim|^|^oQ <f>vix) = 1. Let be a complex- 
valued function on the unit sphere S^, with Jg2 = 1- We use the same 
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symbol for the function on taking values u{x/\x\). For tp and ^ as above, 
consider the expression 

A= dxu{x- yj)VC{x) ■ V(j)v{x - yj) 
JBh 

+ i / v{d{x, yj))ip{x)*(j)y{x - yj)i'{x - yj) . (A.4) 
•IBn 

By using the Cauchy-Schwarz inequality, we can obtain the upper bound 



H'</ dx[\Vi{x)\^ + \v{d{x,y^m{x)\^] (A.5) 
X / dx[\V(l)^{x-yj)\^ + \v{d{x,yj))\(j)^{x-yj)\^]\v{x-yj)\^ . 



I Br 

For an upper bound, we can replace the integration region Br in the second 
integral by M^. Since </>«(x) is a radial function, the angular integration in 
then can be performed by using Jgj \^\'^ = 1- The remaining expression is 
then bounded by a because of /jj3 dx (|V^^(x)p + = 47ra. 

Hence we arrive at 

dx [|V^(x)|2 + \v{d{x,yi)mx)\^] > ^ (A.6) 

Br « 

for any choice of v as above. It remains to derive a lower bound on 
By partial integration, 

/ dx v{x - yj)VC{x) ■ V(f>v{x - yj) 
JBr 

= -/ dxCix)uix - yj)A(t)y{x - yj) 
JBr 

+ dujRCix)i'{x - yj)n-V(f)vix - yj) , (A.7) 

JdBR 

where dujR denotes the surface measure of the boundary of Br, denoted by 
OBr, and n is the outward normal unit vector. Here we used the fact that 
Vz^(x) • V0^(x) = 0. Now, by definition of h{x), ^{x) = il){x) - {2iT)-^/'^h * 
'i/'(x), where * denotes convolution, i.e., h*'ip{x) = Jj^dy h(x — y)'ip{y). Using 
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the zero-energy scattering equation (A. 3) for 0„, we thus see that 



A= duR ip*{x)u{x - yj)n ■ V(f>v{x - yj) 

IdBR 

(27r)-2/2 / <hjR{h*xl)y{x)u{x-yj)n-V(t)v{x-yj) 
JdBn 



+ {2Tr)-^/^ [ dx{h*i^)*{x)u{x-yj)/^(t>y{x-yj). (A.8) 

Jbr 

The last two terms on the right side of (A.8) can be written as 

(27r)-3/2 j [ I df,{y)hiy - x)) 

•JA UBr 



(A.9) 



where d/x is a (non-positive) measure supported in Br. Explicitly, diJ,{x) = 
v{x — yj)A(py{x — yj)dx — n • V(pv{x — yj)v{x — yj)dLOR, the second part 
being supported on the boundary OBr. Note that /^^ dfi = 0, and also 
/g^ c?|/x| < 2a /g2 < 2aV4K (by Schwarz's inequality). Hence 



d/J'iy) h{y - x) 



Br 



< 2aV^7rfR{x -yj) , 



(A.IO) 



with fR defined in (2.6.2). The expression (A.9) is thus bounded from below 

by 

{A.9) > -{2'K)-^lHa^ [ dx\i;{x)\fR{x-yj) 

J A 



>-a[ I dx\ip{x)\ WR{x-yj) 



1/2 



(A.ll) 



Here, we used Schwarz's inequality as well as the definition of wr (2.6.3) in 
the last step. Note that this last expression is independent of ly. 

The only place where ly still enters is the first term on the right side of 
(A.8). By construction, depends only on the direction of the line originat- 
ing from yj, which hits the boundary of Br at a distance not greater than 
R. We distinguish two cases. First, assume that the line hits the boundary 
at a distance R. In this case, we choose ly to be equal to the value of 
at this boundary point. Secondly, if the length of the line is strictly less 
than R, we then choose u to be zero. Of course we also have to normalize 
v appropriately. The integrals are thus only over the part of the boundary 
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of Bj which is at a distance R from yj. Let us denote this part of OBr by 
OBr, assuming for the moment that it is not empty. We then have 



dujRip {x)u{x-yj)n-y^y{x-yj) = R ^ . 

_ (A.12) 
We note that n ■ V(f)v{x — yj) = a/R^ on dBji. We thus obtain from (A. 8)- 
(A.12) 



A>^ dLURmx)fy^ -a(^J^dx\i;{x)\^WRix-yj)y\ (A.13) 



1/2 / /• \ 1/2 

diVR\4^{x)\'' ] 

With the aid of Schwarz's inequality, we see that, for any e > 0, 



\Af>^{l-e) (duR\i;{xt-- [ dx\4^{x)\''wRix-yj). (A.14) 

-K JdBR £ Ja 



At this point wc can also relax the condition that OBr be non-empty; in 
case it is empty, (A.14) holds trivially. 

In combination with (A. 6), (A.14) proves the desired result (A. 2) in 
the special case when C/(|a;|) is a radial (5-function sitting at a radius R, i.e., 
C/(|x|) = R~'^S{\x\ —R). The case of a general potential C/(j.x|) follows simply 
by integrating this result (i.e., Ineq. (A. 2) for this special C/(|x|)) against 
U{R)R'^dR, noting that J dRU{R)R^ < 1 and that wr{x) is pointwise 
monotone increasing in R. 
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